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We calculate the cross sections of electromagnetic productions of muon and tauon pair produc¬ 
tions from the ultra-relativistic heavy ion collisions. Since the Compton wavelengths of muon and 
tauon are comparable to the radius of the colliding ion^ nuclear form factors play important roles 
for calculating the cross sections. Recent measurement indicates that the neutrons are differently 
distributed from the protons therefore this affects the cross section of the heavy lepton pair pro¬ 
duction. In order to see the effects of the neutron distributions in the nucleus, we used analytical 
expression of the Fourier transforms of the Wood-Saxon distribution. Cross section calculations 
show that Wood-Saxon distribution function is more sensitive to the parameter R compare to the 
parameter a. 
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I. INTRODUCTION 

Colliding beams of heavy ions at energies per nucleon in the range 10^ — 10^ GeV, corresponding to experiments 
RHIC at Brookhaven, and LHC at CERN, suggested that intense electromagnetic fields of these ions can produce 
relatively large fluxes of muons and tauons. Since mass of the electron ( 0.511 MeV) is much smaller than muon 
( 105.66 MeV) and tauon (1784 MeV), and the Compton wavelength of the electron (386 fm) is much larger than 
muon (1.86 fm) and tauon (0.11 fm), the cross sections of producing electron pairs are much larger than the heavy 
leptons. In addition to this, since the Compton wavelengths of muons and tauons are smaller than the radius of the 
colliding heavy ions (Au, Pb), nucleon 0 and nucleus form factors are not negligible. Therefore realistic charge form 
factors play important roles for calculating the cross sections of the heavy lepton pair productions. 

In this work, we have used Wood-Saxon nucleus form factor which is also widely used in literature . Instead of 
the monopole approximation, we have an analytical expression of the Wood-Saxon form factor in our equations. With 
the help of this analytical equation, we are able to account for the neutron skin and its effects to our calculations. 
The determination of the sizes and shapes of atomic nuclei is one of the important problems in nuclear physics. 
The rms radius of the charge distribution in nuclei is well known with the high degree of accuracy achieved by the 
elastic electron-nucleus and muon-nucleus scattering experiments. The experimental uncertainties for the nuclear 
charge radius are about 1% for several nuclei Q- However, the neutron distribution and its rms radius in nuclei 
is not precisely determined theoretically and experimentally. Generally, the neutron skin thickness, defined as the 
neutron-proton rms radius difference in the atomic nucleus, 

AR„p = Rn-Rp =< >1/^ -<r^ >y^ (1) 

where the measured charge radius are < 5.43 fm for the gold nucleus and < 5.50 fm for the lead 

nucleus Q. By using the parity-violating electron scattering (PVES), the recent experiment [l| at JLab gives us the 
value Rn as 5.78lQ }g. The difference between the radii of the neutron and proton distributions is the neutron skin 
and it is equal to Rskin = Rn — Rp = 0.33)')Q'itg fm for the lead nucleus with large total error. Therefore the difference 
between the distributions of protons and neutrons in the nucleus affects the cross section of the heavy lepton pair 
production. 

To explore a surface thickness of the neutron density in a heavy nucleus, we model the neutron density with a Wood 
Saxon form: 


Pn{r) = po/[l + exp{r - RQ)/an] ( 2 ) 

where i?o and a„ are the parameters that determine the shape of the neutron distribution. Although the same 
parameters are precisely determined for the proton distribution in the heavy nucleus, there are large uncertainties for 
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FIG. 1: (a) Direct and (b) crossed Feynman diagrams for the heavy lepton pair production in a relativistic heavy ion collision. 

the neutron distribution. In this work, we study the effect of these parameters to the heavy lepton pair production 
cross sections. 

Although the electron-positron pair production from electromagnetic fields goes back to 1930s ii ,there are 
numbers of published literature in the 1990s and 2000s [TMl. The main motivation at these period was the RHIC. 
Dilepton production from the central collisions of heavy ions is very important to understand the hadronic interactions. 
However electromagnetic production of dileptons from the peripheral collisions of heavy ions can shield those dileptons 
produced from the hadronic interactions. Therefore it is highly important to understand various properties of the 
electromagnetically produced dileptons [la¬ 
in our previous works [uJ, Il8|, ll9| , we have calculated electron-positron pair production cross sections by using the 
second order Feynman diagrams. In this work, we have modified the equations by including the electromagnetic form 
factors of nuclei in the momentum space. We have used Fermi type charge form factor. We also consider the coherent 
production of heavy lepton pairs in two photon limit for the symmetric collision of two gold nuclei at RHIC and two 
lead nuclei at LHC energies. In Fig. [U the wiggly lines are the virtual photons that make the electromagnetic field of 
the colliding ions. Here we start with the lowest non-vanishing Feynman diagram which is called two-photon diagram. 
We take the classical limit of the motion of the ions which yields the external field model. Since the parameter 
Za ~ 0.6 is not small. Coulomb and unitarity corrections should be included to obtain exact cross sections. This 
calculations has been done in and we will also mention briefly in the results section. 

In the Formalism section, we briefly explain two photon pair production mechanism. In Calculations and Tables 
section, we study the effects of the form factor that we use in our calculations and we also calculate a variety of cross 
sections for the single-pair process. Finally on the last section, we summarize all the work that we have done in this 
manuscript. 


II. FORMALISM 

In this work we will investigate the two-photon mechanism for heavy lepton pair production. In Fig. [2l the 
electromagnetic interaction of the colliding ions is shown. Since all t he p rotons in the nucleus act coherently, there 
is a very strong electromagnetic field appears around the heavy ions [l0| . Lepton pairs are produced from this field 
which is a cloud of virtual photons. Although we consider peripheral collisions of the heavy-ions where the impact 
parameter 6 > -I-i ?2 is greater than the sum of the radii of the colliding nuclei, our cross section expressions include 
all impact parameter region from zero to infinity. We consider the electromagnetic production of lepton pairs from 
the collisions of heavy ions with charge Z 


Z\ -|- Zi Z\ -\- Zi -f 


where can be electrons, muons or tauons. The maximum magnitude of the electric field of the ions are 


E„ 


Zej 


(3) 

(4) 
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FIG. 2: Collision of two heavy ions at ultra-relativistic velocities. Lorentz contracted strong EM fields are produced around 
the ions. 


where 7 is the Lorentz factor and b is the impact parameter and Ze is the charge of the fully stripped heavy ion. 
For the impact parameters around b ~ 7.5/m which is very closed to nucleus but still peripheral collisions, maximum 
strength of the electric held at the surface E^ax 2 x lO^^V/m for the RHIC and E^ax ^ 7.14 x lO^^R/m for LHC 
collisions. We can estimate the minimum electric held that can create a lepton pair of mass 2mc^ as 


E 


min 


2 4 

me 


ehc 


(5) 


where the distance that we used in the above equation is twice the Compton wavelength of the particular leDton[lll|. 
To create e+e”, and r+r” pairs, the minimum electric held needed are 1.32 x 10^®H/m, 5.66 x lO^^H/m and 

1.62 X respectively. As we see that even at LHC energies, maximum electric held produced from the lead 

ions is below the critical electric held to create tauon pairs. Even for the impact parameters around b ~ 7.5/m —30/m 
regions, maximum electric held is smaller than the electric held that produce tauon pairs at RHIC and LHC energies. 
However, in QED there is still a probability to produce heavy leptons even at these energies. 

The source currents for producing lepton pairs are the Lorentz-boosted charges of the colliding heavy ions. Electro¬ 
magnetic helds in this collisions are very strong that they can pull many lepton pairs from the vacuum. Relativistic 
heavy-ion collider (RHIC) at Brookhaven can accelerate the fully stripped gold nuclei at 100 GeV/nucleon energies 
and large hadron collider (LHC) at CERN can accelerate the fully stripped lead nuclei up to 3400 GeV/nucleon 
energies. Therefore electromagnetic helds at this collisions last for a very short time At = 6 /( 70 /?) and the Fourier 
frequencies could be up to tOmax 'ycPjb [HI, [1^. We can write the Lagrangian density as sum of the non-interacting 
fermion Lagrangian Cn and the Lagrangian density for the coupling of the classical electromagnetic potential to the 
lepton field, mil 


C{x) = Co{x) + Cint{x) 

= 4'(a:)[7;iii9^ — m]4'(a;) — 4'(a:)7^4'(a:)4''(a:) (6) 

where = {Ai,Ai) is the four-vector potential of the nucleus 1 and the scalar component of the potential can be 
written as 

-iq±-b/2 

A° = -8TT^Z-f^S{qo+ Pq;,) GE{q^)fz{q^) (7) 

{qi + Tqi) 

and non-zero component of the vector potential Af is in the longitudinal direction and it can be obtained as 

Al = -PAI ( 8 ) 

where the transverse components of the vector potentials are A\ = A\ = 0. The potentials from nucleus 2 can be 
written by the substitutions, b ^ —b and /? —>■ —/?. In the above equation fz{q^) is the form factor of a nucleus 
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FIG. 3: The electromagnetic form factor for gold. The solid line is the analytic form of Woods-Saxon distribution and the star 
line is the numerical calculation of the same equation. 


which gives the momentum distribution of a proton in the nucleus, and Gsiq^) is the form factor of a proton which 
represents electric distribution of the proton[^[^. As we mention above, we have used two parameter Fermi (2pF) 
function (or Wood-Saxon function) for the charge distribution of the protons: 

= l + eM{r-K)M 

where R is the radius of the nucleus and a is the skin depth or diffuseness parameter. These parameters are obtained 
by fits to electron scattering data and po is written by the normalization condition. For symmetric nuclei, the 
nuclear density for a nucleus that has mass number A, a distance r from its center is modeled in literature with a 
Woods-Saxon distribution as in Eq. ([S]) where po = frn~^ for Au nucleus and po = for Pb nucleus. 

The radii of the gold and lead nucleus are equal to Rau = 6.38 fm and Rpb = 6.62 fm respectively. However, recent 
measurement [I| indicates that the neutrons are differently distributed from the protons so that the parameters R and 
a must be modified. Therefore this effects the cross section of the heavy lepton pair production [^. In order to see 
the effects of the neutron distributions in the nucleus, we need to have analytical expression of the Fourier transforms 
of the Wood-Saxon distribution [2^: 


fz{q^) 



— p{r)Sin{qr)dr 

q 


4 TT 
qr ' 


p{r)Sin{qr)r^ dr 


An'^poa^ 

{qaYSinh‘^{'Kqa) 


OO 

[tt^o Cosh{TTqa)Sin{qR) — qRCos{qR)Sinh^nqa)] + Sirpoa 

n—1 


j^^-nRja 


The parameters R and a are explicitly shown in this analytical equation and it is plotted in Fig. [3] together with 
the numerical result from Fourier transformation of a Woods-Saxon distribution. The agreement is excellent. Since 
the contribution of the last term in above equation is very small compare to the first term, it is safely neglected. By 
changing the parameters R and a we can estimate the sensitivity of the nucleon distributions in the nucleus. 

On the other hand, for the electric distribution function of the proton, we can write the dipole form factor for the 
proton as 

where the value of A^ is 0.71 We have observed that the effects of this distribution is negligible. Finally, 

the total electromagnetic potential is sum of the nucleus 1 and 2 


A'^ = A^ + A^. 


( 12 ) 
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To calculate the total cross section for the heavy lepton pair production, we first represent the direct term, 


= *EE/ i 


du: \Vi>{E^,-^ - UJ) 




(- 


p S 


duj 


- w) 


( 13 ) 


"■hi;/ ZI - £<+>) 

p g J—oc J—oo 

X / d^r Ab{r ; - to) ^ (14) 

J—OO {^p — ‘^) 

Thus, by combining both integrals we obtain for the direct free-free pair production amplitude the explicit expression: 

rflp f dijj , .("l+ii „,b c_2 2 - d"*"’ + - P^)) 


\ <? / J (27r)3 7 27r (fcz - + 7 (kj_ - P±)^ 

S Up 

6(^E^-) -bo- /3(p, - q,)) 1(1 - /3a.)l ui;^) (uip 1(1 + Pa,)\ ui^) 


X 


{pz - +7^(PJ- - Q-l)^ 


ZT-C^) 

rjp — Cc^ 


,(15) 


(s') 

and where Ep is the energy of intermediate state. In these equations k represents the the momentum of electron 
and q represents the the momentum of positron. The vector p describes the momentum of the intermediate (electron 
and positron) states in the field of the ion and can be decomposed into its transverse and parallel part, p = p^ +Pz, 
relative to the motion of the ions. The momenta that is parallel to the heavy ions and the frequency in Eq. m are 
fixed by momentum conservation; 


Pz = 


4"^ - Ei-^^ + I3{k, + q,) 

2P 


(16a) 


_ 4+^ + eE - /3(fc. - 9.) 

2 

For the fixed momentum and the spin states, the transition matrix element can be written as; 


(16b) 


(*7IS.,!*'-') = T / 


cPp± 

(27r)2 




E(ki - Pi : uja) F{p± - Qi : Wb) 7fe5(pi : +13 ), (17) 


uja and oJb are the frequencies associated with the fields of ions a and b, respectively. 


UJa = 


+ /3(9. - kz) 


Wb = 


Ei~^ - 4^^ - I3iqz - K) 


the function E is the scalar part of the field from each heavy ion, 

47rZ7^/3^ 


for the frequency Wq, and as 


E(ki - Pi : Wa) = 


^(Pj. - q± : Wfe) = 


(a; 2 + 72 ^ 2 (kj__p^) 2 ) 1 

47rZe7^/3^ 

R+ 4 ¥^(p 7 ^TE) ’ 


(18a) 

(18b) 

(19a) 

(19b) 
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for the frequency ojb, respectively. The function 7” explicitly depends on the velocity of the heavy ions (3, on the 
transverse momentum pj_, and on the states k,q; 


rkq{p± : +/3) = 5ZX! 

S (Jp 



X 1(1 1(1 . 

Having the amplitudes for the direct and crossed diagram, we are now prepared to 
the generation of a free-free pair in collisions of two heavy ions 

•' k>0 q<0 


( 20 ) 

write down the cross section for 

( 21 ) 


where S = Sab + Sba denotes the sum of the direct and crossed. Making use of all the simplifications from above, 
these cross sections can be expressed as: 


with 


and 


k>0 q<0 

“Hi:/ (.4W(t,,;px)+.4l-)(C«pU))b (22) 

<yk <7q ^ ' 

A^^\k,q-,p±) = F(k± - p± : coa) F{p± - q± : uJb)Tkq{p± ■■+13), (23a) 

g;p_L) = F{k± - p± : uJb) F{p± - c[± : u}a)Tkq{p± ■ (23b) 


being some proper products of the transition amplitudes and scalar parts of the fields as associated with ions a and b. 


III. CALCULATIONS AND TABLES 

Table U displays the cross sections of electromagnetic productions of electron, muon and tauon pairs from the ultra- 
relativistic gold-gold collisions at RHIC energies and lead-lead collisions at LHC energies with no form factor and 
with Wood-Saxon form factor distribution. In this calculation, we simply used the analytical expression of the Fourier 
transformation of this form factor. We have also include the electric distribution function of the proton (Eq. 9) in 
our calculation, and we have seen that it has negligible effect for the production of any lepton pairs. 

When we compare electron pair production cross sections with and without form factor, it is clearly seen that the 
difference among the results are negligibly small. Since the Compton wavelength of the electron is much larger than 
the radius of the colliding heavy ions, form factors of the proton and nucleus almost have no effect for the production 
of electron-positron pairs. 

However, when we compare muon and tauon pair production cross sections with and without form factors, it is 
clearly seen that the differences are large. The cross section results for muon and tauon pair production with form 
factors are much smaller than without form factors. Since, the Compton wavelengths of muons and tauons are smaller 
than the radius of the colliding heavy ions, the effects of the form factors become dominant for pair production. At 
RHIC and LHC energies, muon pair production is reduced by about 3 and 2 factors, respectively. On the other hand, 
at RHIC and LHC energies, tauon pair production is reduced by about 100 and 5 factors, respectively. 

Figs, m [5] and m displays the differential cross sections as a function of the energy po, the longitudinal momentum 
Pz, the transverse momentum p± and rapidity y of the produced electron, muon and tauon for with and without 
form factors at RHIC and LHC energies. When we compare these figures, in Fig H] for the electron pairs, it is clearly 
seen that form factor has almost no effect to the cross section. With and without form factors, differential cross 
sections have the same behavior and almost identical values. The differential cross sections as a function of transverse 
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TABLE I: The cross sections of electromagnetic productions of electron, muon and tauon pair productions from the ultra- 
relativistic gold-gold collisions at RHIC energies and lead-lead collisions at LHC energies. means form factor is not 

included and means that form factor is included in the calculations. 


fNFF 

(barn) 


fWS 

cr^ 


(barn) 


RHIC 

3.63 X 10^ 

3.62 X 10' 

LHC 

2.44 X 10® 

2.43 X 10' 

RHIC 

1 

o 

t-H 

X 

CO 

2.03 X 10’ 

LHC 

5.2 

2.7 

RHIC 

4.7 X lO""* 

3.9 X 10“' 

LHC 

7.7 X 10"® 

1.6 X 10“ 


TABLE If: The percentage of cross section contributions to the total cross section between the impact parameter regions 
0 < 6 < 15 — 20 — 25 — 30/m and 15 — 20 — 25 — 30/m < b < oo. The cross section calculations are taken from the Table I with 
form factor. 



e 

b < 15/m 

'6+ 

b > 15/m 

b < 15/m 

+ 

b > 15/m 

T 

b < 15/m 

T + 

b > 15/m 

RHIC 

LHC 

0.04% 

0.04% 

99.96% 

99.96% 

13.2% 

3.65% 

86.8% 

96.35% 

49.8% 

25% 

50.2% 

75% 


b < 20 fm 

b > 20fm 

b < 20fm 

b > 20fm 

b < 20fm 

b > 20fm 

RHIC 

LHC 

0.07% 

0.07% 

99.93% 

99.93% 

20.44% 

6.24% 

79.56% 

93.76% 

60.11% 

35.23% 

39.89% 

64.77% 


b < 25/m 

b > 25/m 

b < 25/m 

b > 25fm 

b < 25fm 

b > 25/m 

RHIC 

LHC 

0.11% 

0.11% 

99.89% 

99.89% 

27.6% 

9.27% 

72.4% 

90.73% 

67.13% 

43.77% 

32.87% 

56.23% 


b < 30/m 

b > 30fm 

b < 30fm 

b > 30fm 

b < 30/m 

b > 30fm 

RHIC 

LHC 

0.16% 

0.16% 

99.84% 

99.84% 

34.15% 

12.6% 

65.85% 

87.4% 

72.15% 

50.7% 

27.85% 

49.3% 


momentum decreases rapidly compare to the longitudinal momentum and energy. Differential cross sections of as 
a function longitudinal momentum and energy have almost identical values. This also shows that produced leptons 
carry energy in the longitudinal direction where the heavy-ions move. 

On the other hand, In FiglS] and [5] the same differential cross sections of the produced muon and tauon pairs at 
RHIC and LHC energies with and without form factors are shown. When we compare these figures, it is clear that the 
longitudinal momentum of the produced muon and tauon are much higher than the transverse momentum. Therefore, 
the energy E = are carried by the longitudinal momentum of the produced heavy leptons. It is also clear 

that, the differential cross sections as a function of energy and longitudinal momentum are very close to each other 
at RHIC and LHC energies. The rapidity distribution of muon pair production, form factor reduce the cross sections 
by about ten times at the RHIC energies, and reduce the cross section by about five times at the LHC energies. For 
the tauon pair production, the reduction of the cross section is about two order of magnitude for the RHIC energies, 
and is about one order of magnitude for the LHC energies. All these calculations show that distributions of the heavy 
leptons are mainly confined to the central rapidity region. 
















TABLE III: The ratio of the cross sections (Acr/a) of the produced muon and tauon pairs. At the first column, the parameter 
R is kept constant and the parameter a is changed between 0.4/m < a < 0.6fm (Aa = 0.2/m). At the second column, the 
parameter a = 0.55fm is kept constant and the radius R is changed between 6.1 fm < R < 7.0/m. 


Rau 

= 6.38fm,Rpi, = 6.62/m Aa = 0.2/m 
Aa j a 

a — 0.55/m, AR — 0.9/m 
Aa f a 

p-p+ RHIC 

4% 

11 % 

LHC 

2 % 

5 % 

r'r+ RHIC 

14 % 

35 % 

LHC 

4 % 

10 % 






FIG. 4: The differential cross section as function of a) energy (po), b) longitudinal momentum (pz), c) transverse momentum 
(p±) and d) rapidity (y) of the produced electron at RHIC and LHC energies with and without form factors. 


Fig. [7] displays muon and tauon total cross sections for Au + Au and Pb + Pb collisions as a function of the Lorentz 
factor 7 of the beams for the Wood-Saxon form factors. We also compare this results with the no form factor. For the 
muon pair production, the total cross sections are reduced about 10 factors for the low energies, about 3 factors for 
the RHIC energies, and about 2 factors for the LHC energies. For the tauon pair production, the total cross sections 
are reduced about 10® — 10® factors for the low energies, about 100 factors for the RHIC energies, and finally about 
5 factors for the LHC energies. 

Fig. [5] shows the total cross section of heavy lepton pair production as a function of the charge of the colliding 
nuclei Z. Obviously nuclear form factor reduces the total cross sections. At RHIC energies, the form factor reduces 
the tauon pair production cross section about one or two order of magnitude, however for the muon pair production 
the reduction is about two or three factors. On the other hand, at LHC energies, form factors reduce the muon pair 
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(c) 


(d) 


FIG. 5: The differential cross section as function of a) energy (po), b) longitudinal momentum (pz), c) transverse momentum 
(p_l) and d) rapidity (y) of the produced muon at RHIC and LHC energies with and without form factors. 


production cross section about one or two factors, for the tauon production cross section this reduction is slightly less 
than one order of magnitude. 

Fig. [9] displays the impact parameter dependence cross section for muon and tauon pairs production at RHIC and 
LHC energies with and without Wood-Saxon distribution. In this calculations, the range of the impact parameter b 
is between 0 and infinity. Since the Compton wavelengths of muon is 1. 86 /m and tauon is O.II/m which are smaller 
than the radius of the nucleus, the heavy leptons are produced mainly near the radius of the heavy ions. Maximum 
production of muon pairs occurs at RHIC energies about impact parameters of 20 fm region, and at LHC energies 
in the impact parameter of 40 - 80 fm region. At about 200 fm region, the impact parameter dependence cross 
section decreases at about 10 factors at RHIC energies, and 2 factors at LHC energies. On the other hand, maximum 
production of tauon pairs occurs about impact parameters of 7 fm region at RHIC energies, and about 10 fm region 
at LHC energies. At about 200 fm, the impact parameter dependence cross section decreases about 100 factors at 
RHIC energies and 50 factors at LHC energies. 

Table m represents the percentage of cross section contributions to the total cross section between the impact 
parameter regions b < Ri + R 2 and b > Ri + R 2 where i?i + i ?2 ~ 15/m for the heavy ions. The cross section 
calculations are taken from the Table I with form factor. For the electron-positron pair production cross section, pairs 
are overwhelmingly produced for the region of 6 > i?i -I- i? 2 - For the muon and tauon pair production, b < Ri + R 2 
region has important contribution to the total cross sections. When the energy of the colliding ions increases, the 
contributions from the b < R 1 +R 2 region decreases. However, when the mass of the lepton increases, the contribution 
from the b < Ri + R 2 region increases. We have also repeated the calculations for the impact parameters b for the 
regions 20, 25 and 30 fm. This calculations clearly indicate that the details of the form factors, such as the distributions 
of the nucleons in the nucleus, are important for the calculations of the cross sections. 

In all above calculations, we have also used proton form factor GE{q^)- Our results show that the contribution of 
this proton form factor is negligible. 

The description of matter distribution in nuclei is still under investigation in nuclear physics and precise knowledge 
of this distribution can help us to understand the structure of the nucleus. Although the charge distribution has 
been measured with high accuracy from the electron scattering data so that the charge radii are usually known with 
uncertainties lower than 1% [1,2]. However, the neutron distribution in the nucleus is not accurately measured as 






















10 






FIG. 6: The differential cross section as function of a) energy (po), b) longitudinal momentum (pz), c) transverse momentum 
(p_l) and d) rapidity (t/) of the produced tauon at RHIC and LHC energies with and without form factors. 


protons. 

We have also investigated the neutron distribution of the nucleus. Wood Saxon equation gives us the proton 
charge distribution and neutron density distribution of the nucleus. In Eq. 8, the parameter i? is the radius of the 
nucleus and a is the skin depth and these parameters are obtained from the electron scattering data [^ . Recent 
experiments show that neutrons are distributed differently from the protons and this effects the cross sections of the 
heavy lepton pair productions. Neutrons are electrically neutral, so do not contribute to the electric form factor. 
Neutrons can contribute to the hadronic radius. In order to see the effects of the neutron distribution, we have 
changed the parameters R and a in the Wood-Saxon form factor. First we keep the R constant and change the 
parameter a between 0.3 /to < a < 0.8/to. Although the parameter a has a value of about 0.55 /to, by changing it for 
the extreme ranges, we can see how it effects the cross sections. Next we keep the parameter a = 0.55/to constant and 
the change the radius R between 5.8 /to < R < 7.2/to. In Figs. ITUlandfm we have plotted these cross sections for 
muon and tauon productions at RHIC (gold nucleus) and at LHC (lead nucleus). When the parameter a is changed 
from 0.4 /to to 0.6fm value, the cross section (Acr/tr) of the produced muon pairs decreased by about 4% at RHIC 
and 2% percent at LHC. Similarly, tauon production decreased by about 4% at LHC and 14% at RHIC. On the other 
hand, when we change the radius R from from 6.1 fm to 7.0 /to, the cross sections of the muon pairs decreased by 
about 5% at LHC and 11% at RHIC. The similar calculation shows that tauon production decreased by about 10% at 
LHC and 35% at RHIC. These calculations are tabulated in Table HI. From this calculations, we can conclude that 
Wood-Saxon distribution function is more sensitive to the parameter R compare to the parameter a. It also shows 
that, cross sections are very sensitive to the parameter R at RHIC energies compare to the LHC energies. The proton 
distribution in the nucleus is well known. If the neutron surface is shifted outwards the proton surface this means 
that the skin depths of neutrons and protons are nearly equal a„ ~ Op and it is called bulk or skin type of neutron 
skin. On the other hand, if the diffuseness of the neutron surface is greater than the proton then ^ Rp and it is 
called surface or halo type of the neutron skin [30l |. 

We have compared our results with the works done in Ref.[3l|. Both cross section results of the RHIC for 
production nearly agree each other, however cross section in our work approximately is 3 times larger than 

their results. Our LHC cross section results of IJ .~and r“T+ productions are slightly higher than their results. 
Authors in Ref. P they used Caussian distribution form factor and the equivalent photon approximation to calculate 
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(C) (d) 


FIG. 7: Muon and tauon cross sections for Au + Au and Pb + Pb collisions as function of the 7 with and without form factors. 


the heavy lepton pair production cross section. 

When we compare our cross section results with Ref. [ 2 ^, both calculations for fJ,~ production at RHIC energies 
are almost similar with and without form factor At RHIC energies, r“r+ cross section with no form factor in our 
work approximately is 5 times lower than in Ref. [ 2 ^, however the results with form factor is close to each other. 

Our cross section results of rates at RHIC and LHC are about slightly lower than in Ref. [ 2 ^ that is used 

Born calculations at RHIC and LHC energies. Since the parameter Za « 0.6 is not small, Coulomb and unitarity 
corrections 


^total — ^two—photon ^Coulomb “I ^unitary 

should be included in the calculations to obtain exact cross sections. In Ref. [ 2 ^ authors show that for electron positron 
pair production, the Coulomb corrections to the Born cross section are large, whereas the unitary corrections are small. 
On the other hand, for muon pair production the Coulomb corrections to the Born cross section (corresponding 
to multiphoton exchange of the produced muon pairs with the nuclei) are small, whereas the unitary corrections 
(corresponding to the exchange of light-by-light blocks between nuclei) are large. Authors find that for Pb + Pb 
collisions at LHC there is about —50% unitarity corrections to the Born cross sections. 


IV. RESULTS 

In this work, we investigate the contribution of the Wood-Saxon nuclear form factors on the heavy lepton pair 
production. We study the pair productions for the RHIC and LHC energies for the Au and Pb nucleus. Even though 
electron pair production is not strongly effected by the form factors, the calculations show that nuclear form factors 
play important role for the muon and tauon pair production. At RHIC and LHC energies, muon pair production is 
reduced by about 3 and 2 factors by the Au and Pb form factors, respectively. On the other hand, at RHIC and LHC 
energies, tauon pair production is reduced by about 100 and 5 factors by the Au and Pb form factors, respectively. 
These calculations show that tauon pair production is very sensitive to the form factors at the RHIC energies. We 
can conclude that, muon pair production in peripheral heavy-ion collisions is coherent greater than the dimensions 
of the nucleus. However, although for the LHC energies the tauon pair production is coherent over distances of the 
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(a) (b) 




(c) (d) 

FIG. 8: Muon and tauon cross sections for RHIC and LHC energies as functions of the nuclear charge Z with and without form 
factors. 




(a) (b) 

FIG. 9: Impact parameter dependence cross section for muon and tauon pairs at RHIC and LHC energies with Wood-Saxon 
distribution and no form factor. 


dimension of the nucleus, for the RHIC energies this coherence is break down. Since the mass of the tauon is very 
large when compared to muon, incoherent pair production of tauon is dominant at the RHIC energies. When the 
energy is increased, the coherence pair production is restored. We expect that this result is also valid for the other 
heavy particles. 

We have also investigated the effect of the neutron distributions in the nucleus to the cross section of the heavy 
lepton pair production. Neutron skin can be produced in two ways: First type is as a result of the shift of the neutron 
surface in relation to the proton one (the bulk or ’’skin” type a„ ^ Up), and the second type is due to the differences 
in the surface thickness of proton and neutron matter in a nucleus (the surface or ’’halo” type ~ i?p). In general, 
both types of mechanisms contribute to the creation of neutron skin. We have changed the parameters a and H in 
the Wood-Saxon form factor, and our calculations show that Wood-Saxon distribution function is more sensitive to 
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a{fm) R{fm) 

(C) (d) 

FIG. 10: The cross section of produced muon a) as a function of skin thickness a at RHIC energy when R = 6.38/m b) as a 
function of R at RHIC energy when skin thickness a = 0.55/m. The cross section of produced tauon c) as a function of skin 
thickness a at RHIC energy when R = 6.38/m d) as a function of R at RHIC energy when skin thickness a = 0.55/m. 


the parameter R compare to the parameter a. 
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